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Variational weak-coupling perturbation theory yields converging approx- 
imations, uniformly in the coupling strength. This allows us to calculate 
directly the coefficients of strong- coupling expansions. For the anharmonic 
oscillator we explain the physical origin of the empirically observed conver- 
gence behavior which is exponentially fast with superimposed oscillations. 

1 Introduction 

An important problem of perturbation theory is the calculation of physically mean- 
ingful numbers from expansions which are usually divergent asymptotic series with 
coefficients growing oc k\ in high orders k. For small expansion parameters g a 
direct evaluation of the series truncated at a finite order k « 1/g can yield a rea- 
sonably good approximation, but for larger couplings such series become completely 
useless and require some kind of resummation. Well-known examples are field the- 
oretical e-expansions for the computation of critical exponents of phase transitions, 
but also the standard Stark and Zeeman effects in atomic physics lead to divergent 
perturbation expansions. 

The paradigm for studying this problem is the quantum mechanical anhar- 
monic oscillator with a potential V{x) = \ui 2 x 2 + \gx 4 (u 2 ,g > 0). The Rayleigh- 
Schrodinger perturbation theory yields for the ground-state energy a power-series 
expansion 



where the E£> are rational numbers 1/2, 3/4, -21/8, 333/16, -30885/128, ... , 
which can easily be obtained to very high orders from the recursion relations of 
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Bender and Wu [1]. Their large-order behavior is analytically known to exhibit the 
typical factorial growth, 

4 0) = -aM(QM l/2 (-^) k k~ 1/2 k\(i + 0{l/k)). (2) 

Standard resummation methods are Pade or Borel techniques whose accuracy, 
however, decreases rapidly in the strong-coupling limit. In this note we sum- 
marize recent work on a new approach based on variational perturbation theory 
[2, 3]. Our results demonstrate that by this means the divergent series expan- 
sion (1) can be converted into a sequence of exponentially fast converging ap- 
proximations, uniformly in the coupling strength g [4-7]. This allows us to take 
all expressions directly to the strong-coupling limit, yielding a simple scheme for 
calculating the coefficients aij of the convergent strong-coupling series expansion, 
E(°\g) = (g/A) 1 / 3 \a + a^Au 3 / g) 2 ' 3 + a 2 {Au 3 / gfl 3 + . . .} . 



2 Variational Perturbation Theory 

The origin of variational perturbation theory can be traced back to a variational 
principle for the evaluation of quantum partition functions in the path-integral for- 
mulation [3, 8]. While in many applications the accuracy was found to be excellent 
over a wide range of temperatures, slight deviations from exact or simulation results 
at very low temperatures motivated a systematic study of higher-order corrections 
[2, 3]. 

In the zero-temperature limit the calculations simplify and lead to a resumma- 
tion scheme for the energy eigenvalues which can be summarized as follows. First, 
the harmonic term of the potential is split into a new harmonic term with a trial 
frequency f2 and a remainder, uj 2 x 2 = Vl 2 x 2 + (oo 2 — il 2 )x 2 , and the potential is 
rewritten as V(x) = \VL 2 x 2 + \g{— 2ax 2 /Q + x A ), where a = Q(Q 2 — uj 2 )/g. One 
then performs a perturbation expansion in powers of g = g/Q 3 at a fixed a, 

N 

^WHE^VM^f, (3) 

fc=0 

where E$ = E^/Q is the dimensionless reduced energy. The new expansion 
coefficients are easily found by inserting u = ^JQ 2 — ga/Q = — ga in (1) 
and reexpanding in powers of g, 

4 0) W = E< ) ( (1 ^ )/2 )(-4^. (4) 

The truncated power series Wn(9, ^) = (g,cr) is certainly independent of 

Q in the limit iV — > oo. At any finite order, however, it does depend on Q, the 
approximation having its fastest speed of convergence where it depends least on Q, 
i.e., at points where 3W N /dVt = 0. If we denote the order-dependent optimal value 
of Vt by Qjv, the quantity Wjv(<7,^jv) lS the new approximation to E^\g). 
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At first sight the extremization condition OWn/OQ = seems to require the 
determination of the roots of a polynomial in of degree 3N, separately for each 
value of g. In Ref. [4] we observed, however, that this task can be greatly simplified. 
While Wn does depend on both g and Q separately, we could prove that the deriva- 
tive can be written as OW^/dVt = (g/4) N P N (a), where Pjv(ct) = —2de$ +1 (a) / da is 
a polynomial of degree N in a. The optimal values of a were found to be well fitted 
by 

a N = cN (l + 6.85/iV 2/3 ) , (5) 

with c = 0.186 047272... determined analytically (cp. Sec. 3). This observation 
simplifies the calculations considerably and shows that the optimal solutions fl^ 
depend only trivially on g through a^ = Q^(Q 2 N — uo 2 )/ g. Since the explicit knowl- 



edge of f2jv is only needed in the final step when going back from E^' to E y N . 



this suggests that the variational resummation scheme can be taken directly to the 
strong-coupling limit. 

To this end we introduce the reduced frequency Cj — u/Q, write the approxima- 
tion as Wn = (g/g) 1 ^ 3 u>n(£I, ^ 2 ), and expand the function WN(g,& 2 ) in powers of 



»/0)2/3£2/ 3- Xhig giveg [5] 



W N = (g/4) 
with the coefficients, 



1/3 



a + «! (4uj 3 / g y /3 + a 2 (4cj 3 /^) 



4/3 



+ . . . 



(6) 



OLr. 



(2ra-l)/3 




-g/^Y- (7) 



If this is evaluated at g — 1/crjv with a^ given in (5), we obtain the exponentially 
fast approach to the exact limit as shown in Fig. 1 for olq. The exponential falloff is 
modulated by oscillations. Our result, a = 0.667986 259 155 777108 270 96, agrees 
to all 23 digits with the most accurate 62-digit value in the literature. The compu- 
tation of the higher-order coefficients a n for n > proceeds similarly and the results 
up to n = 22 are given in Table 1 of Ref. [5]. 



3 Convergence Behavior 

To explain the convergence behavior [6, 7] we recall that the ground-state energy 
(g) satisfies a subtracted dispersion relation which leads to an integral represen- 
tation of the original perturbation coefficients, 

where disc E^ (g) = 2ilmE^ {g — irj) denotes the discontinuity across the left-hand 
cut in the complex (/-plane. For large k, only its g — > 0~ behavior is relevant and a 
semiclassical calculation yields disc E^\g) m — 2iuj(6/n) l / 2 (— 4ui 3 /3g) 1 ^ 2 exp(4u; 3 /3g), 
which in turn implies the large-order behavior (2) of Ef\ 
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Figure 1: L.h.s.: Exponentially fast convergence of the Nth approximants for a to 
the exact value. The dots show Ajv = Kcko)^ — «o|- Rh.s.: Cuts in the complex 
g-plane. The cuts inside the shaded circle happen to be absent due to the converge 
of the strong-coupling expansion for g > g s . 

The reexpanded series (3) is obtained from (1) by the replacement of uo — > 
— erg. In terms of the coupling constant, the above replacement amounts to 
g = g/uj 3 — ► g/(l — ag) 3 / 2 . Using this mapping it is straightforward to show [6] that 
£■(0) = E^/Q satisfies a dispersion relation in the complex g-plane. If C denotes 
the cuts in this plane and discc-E*- ^^) is the discontinuity across these cuts, the 
dispersion integral for the expansion coefficients reads 



In the complex g-plane, the cuts C run along the contours Ci, C\, C 2 , Cg, and C 3 , 
as shown on the r.h.s. of Fig. 1. The first four cuts are the images of the left-hand 
cut in the complex g-plane, and the curve C3 is due to the square root of 1 — ag in 
the mapping from g to g. 

Let us now discuss the contributions of the various cuts to the fcth term Sk- For 
the cut C\ and the empirically observed optimal solutions a N = cN(l + b/N 2 / 3 ), a 
saddle-point approximation shows [6] that this term gives a convergent contribution, 
SW(Ci) oc e -[- 61 °s(-T)+( c s)" 2/3 ] Jvl/3 , only if one chooses c = 0.186 047272... and 
7 = —0.242 964 029 . . .. Inserting the fitted value of b = 6.85 this yields an exponent 
of — 61og(— 7) = 9.7, in rough agreement with the convergence seen in Fig. 1. If 
this was the only contribution the convergence behavior could be changed at will by 
varying the parameter b. For b < 6.85, a slower convergence was indeed observed. 
The convergence cannot be improved, however, by choosing b > 6.85, since the 
optimal convergence is limited by the contributions of the other cuts. 

The cut Ci is still harmless; it contributes a last term S N (Ci) of the negligible 
order e - Nl °s N , The cuts C 2 ,2,3, however, deserve a careful consideration. If they 
would really start at g — 1/a, the leading behavior would be £^(62,2,3) °c °"* \ an d 
therefore Sn(C 2> 2,3) oc (crg) N , which would be in contradiction to the empirically 




(9) 
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observed convergence in the strong-coupling limit. The important point is that the 
cuts in Fig. 1 do not really reach the point ag — 1. There exists a small circle of 
radius Ag > in which (g) has no singularities at all, a consequence of the fact 
that the strong-coupling expansion (6) converges for g > g s . The complex conjugate 
pair of singularities gives a contribution, 

S N (C 2X3 ) w e- Nl/aacos9 cos(N 1/3 asm6), (10) 

with a = 1/ (\g s \c) 2 / 3 . By analyzing the convergence behavior of the strong-coupling 
series we find \g s \ ^0.160 and 8~— 0.467, which implies for the envelope an asymp- 
totic falloff of e ~ 9 - 23Nl/s , and furthermore also explains the oscillations in the data [6] . 

4 Conclusions 

To summarize, we have shown how variational perturbation theory can be used to 
convert the divergent weak-coupling perturbation series of the anharmonic oscillator 
into a sequence of converging approximations for the strong-coupling expansion. By 
making use of dispersion relations and identifying the relevant singularities we are 
able to explain the exponentially fast convergence with superimposed oscillations in 
the strong-coupling limit. 

W.J. thanks the Deutsche Forschungsgemeinschaft for a Heisenberg fellowship. 

References 

[1] CM. Bender and T.T. Wu, Phys. Rev. 184 (1969) 1231; Phys. Rev. D7 (1973) 
1620. 

[2] H. Kleinert, Phys. Lett. A173 (1993) 332; H. Kleinert and H. Meyer, Phys. 
Lett. A184 (1994) 319. 

[3] H. Kleinert, Path Integrals in Quantum Mechanics, Statistics and Polymer 
Physics, 2nd edition (World Scientific, Singapore, 1995). 

[4] W. Janke and H. Kleinert, Phys. Lett. A199 (1995) 287. 

[5] W. Janke and H. Kleinert, Phys. Rev. Lett. 75 (1995) 2787. 

[6] H. Kleinert and W. Janke, Phys. Lett. A206 (1995) 283. 

[7] For related work, see I.R.C Buckley, A. Duncan, and H.F. Jones, Phys. Rev. 
D47 (1993) 2554; A. Duncan and H.F. Jones, Phys. Rev. D47 (1993) 2560; 
CM. Bender, A. Duncan, and H.F. Jones, Phys. Rev. D49 (1994) 4219; C 
Arvanitis, H.F. Jones, and CS. Parker, Phys. Rev. D52 (1995) 3704; R. Guida, 
K. Konishi, and H. Suzuki, Ann. Phys. 241 (1995) 152; 249 (1996) 109. 

[8] R.P Feynman and H. Kleinert, Phys. Rev. A34 (1986) 5080; R. Giachetti and 
V. Tognetti, Phys. Rev. Lett. 55 (1985) 912; Phys. Rev. B33 (1986) 7647. 



5 



